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Abstract 

I discuss the quantities and effects important for the freeze-out and outline 
a formalism for the description of continuous decoupling of particles from 
the fireball. Then I present a calculation of the scattering rates of pious at 
^ , various temperatures and argue that it is important to take continuous particle 

■ decoupling into account when modelling the collision dynamics. 

Hadronic single-particle spectra and correlations carry information about the 
(<~^ , freeze-out state of the fireball which results from the collective expansion of the 

strongly interacting matter. Often, the freeze-out state is being modelled in frame- 
work of the Cooper- Frye ^ mechanism where all particles — regardless their identity 
and/or momentum — are emitted from a single sharp three-dimensional freeze-out 
O , hypersurface. This is the case for most hydrodynamic simulations. A question ap- 

pears: does the assumption of a common sharp three-dimensional freeze-out still 
provide reliable approximation of the real process? 

In case of a sharp freeze-out, the freeze-out hypersurface is usually characterised 
^> , by some prescription. There were some attempts to identify a universal freeze-out 

^ • criterion i-e., a condition which determines the freeze-out hypersurface for 

5^ , heavy-ion collisions of any size and at any ultrarelativistic energy. In case of the 

continuous gradual decoupling there is no hypersurface to be determined and the 
concept of a universal freeze-out criterion is not applicable. I will focus on gradual 
decoupling since the sharp freeze-out can be defined as its limiting case. 

Let me focus on the mechanism of freeze-out and identify the important effects 
and quantities. Freeze-out occurs when scattering ceases. It has been suggested 
that the mean- free path is the relevant quantity to look at Here, densities of 

the individual species are weighted with cross-sections for scattering on them. An 
example: density of nucleons is much more important for pion scattering rate than 
the density of pions, because the ttN cross-section is bigger than the one for tttt 
scattering. The CERES collaboration argued that the "universal" pion mean-free 
path at freeze-out should be something of the order of 1 fm (maybe 2-3 fm) |3 . 

On the first sight, this is rather surprising, because this length is much shorter 
than the size of the system. So far, however, we did not mention the expansion 
and the decrease of the density due to it. In rough terms, freeze-out occurs when 
the rate of the density decrease becomes comparable or larger than the scattering 
rate [H]. The picture is the following: at the moment of scattering the particle has 
some finite mean-free path A and thus would be expected to scatter after some time. 
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However, if the density drops fast enough, after passing the length A our particle 
may find itself in an environment of very low density such that no other scattering 
occurs. 

In [S] we used the formalism of escape probabilities (conceived earlier by other 
authors, e.g. [3111) and focused on the dependence of the escape probability on 
temperature and the momentum of pions. By using this formalism one qualitatively 
goes beyond probing the mean-free path and the rate of density decrease; particles 
now decouple from a finite four-volume. Sharp freeze-out is realised as a limiting 
case of this more general case, when the escape probabilities of all particles change 
from to 1 in a very narrow space-time region. 

The escape probability can be determined as 

V{x,T,p) = exp ^— J df TZ{x + vf , p)^ , (1) 

where TZ is the scattering rate. The opacity integral in this equation is evaluated 
along the expected trajectory of a particle with momentum p and the position char- 
acterised by X and t if it moved straight. The integrated scattering rate gives the 
average number of collisions the particle would suffer on this trajectory. Assump- 
tions about the chemical composition of the medium are included in calculation of 
TZ, while the density decrease rate determines its time dependence and thus the 
value of the resulting integral. 

As an example of the fireball expansion we assumed in 6 a particle with pt — 
in the centre of the fireball which does not move transversely but just waits until 
the surrounding matter decays. Then we only need a prescription for the time 
dependence of TZ for which we take an ansatz: 

7^(r)=7^o(^)^ (2) 

where a > 1 and TZq is the scattering rate at time tq. This corresponds to a fireball 
expanding longitudinally in a boost-invariant manner and transversely with the 
following dependence of the transverse flow rapidity on radial coordinate r: 

a — 1 r 

Vt\r,V ^ (3) 
The escape probability in this toy model then reads 

Vip, = 0,r = 0) = exp f-^) = exp (- ) . (4) 

Obviously, larger TZq leads to smaller escape probability, whereas stronger expansion 
is encoded in larger a and/or r^t and increases the value of "P. 

We calculated the scattering rate of negative pions in thermal hadronic gas as 

n{p)^Y. J d^kp^{k)al{s)\v^-v,\* . (5) 

i 

Here, pi{k) is the density of particles on which the pion scatters, <j'^{s) is the collinear 
cross-section, and — | * is the relative velocity of the pion to the other scattering 
partner in the CMS of the pair. We integrate over the momenta of the scatterers 
and sum over following species as scattering partners: i — tt, N, TV, K, p, A, A. 
For the cross-section we use a parametrisation of hadron scattering via resonances 
(see iniini). For the temperatures we assumed values of 90, 100, and 120 MeV. We 
estimated the chemical potentials in such a way that we reproduced data on pion 
phase-space density at the SPS |10l lll| and RHIC and the ratios of dN/dy of 
different species at midrapidity [T^IT^IT^] . 
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Figure 1: Scattering rate as a function of pion momentum with respect to the 
hadronic medium with T = 100 MeV. Results are obtained for SPS (left) and RHIC 
(right). Contributions to the scattering rate from scattering on pions, nucleons and 
antinucleons are indicated. The two lower panels show the baryonic and mesonic 
relative contributions. 



In Fig. n I compare pion scattering rates in hadronic gases corresponding to 
those at the SPS and at RHIC. The contribution from nucleons and antinucleons to 
the total scattering rate increases little when moving from SPS to the higher energy 
system. A higher phase-space density of pions at RHIC leads to larger pion 
contribution, but the total scattering rate is not dominated by pions and therefore 
does not change much. Quantitatively, about 10% relative contribution is shifted 
from baryons to mesons when going from SPS to RHIC. 

A summary of results is plotted in Fig. 12 The scattering rate drops with temper- 
ature. When we use eq. Q and a realistic estimate of the transverse rapidity gradi- 
ent ^B], reasonable escape probabilities (~ 30-50%) are obtained for T < 100 MeV. 
Note that this was inferred for a particle with = in centre of the fireball. 

One can see that regardless the assumed temperature and chemical potentials 
higher momentum particles always have smaller scattering rate and thus decouple 
easier. This suggests that the freeze-out can be sequential: high momenta first, low 
momenta later. A sequential freeze-out might be an additional cause of the Mt 
dependence of the HBT radii, if high-pt pions decouple earlier and from a smaller 
fireball than the low-pt ones. Such a scenario, however, must be studied in more 
detail. Nevertheless, it represents an interesting alternative to the blast-wave model, 
which seems to have problems in reproducing data on HBT radii |17j . 

Due to the momentum dependence of the scattering rate, collapsing the whole 
decoupling four- volume into a sharp three-dimensional freeze-out hypersurface seems 
an unreliable approximation. 
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Figure 2: The pion scattering rates as functions of momentum with respect to 
the medium, calculated for different temperatures and sets of chemical potentials 
allowed by data. Details can be found in the original paper . 
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